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ABSTRACT 

We analyze differen t possibilities to expla in the wide initial spin period distribution of radio 
pulsars presented bv lNoutsos et al.l(l2013l) . With a population synthesis modeling we demon- 
strate tha t magnetic field deca y can be use d to interpret the differe nce between the recent 
results bv lNoutsos et al.l (1201 3l) and those by lPopov & Turolial(l2012h . where a much younger 
population of NSs associated with supernova remnants with known ages has been studied. In 
particular, an exponential field decay with r^ag = 5 M yrs can produce a "ta il" in the recon- 
structed initial spin period distribution (as obtained by iNoutsos et al.ll2013b up to f o > Is 
starting with a standard gaussian with (Pq) = 0.3 s and crp - 0.1 5 s. Another option to ex- 
plain the difference between initial spin period distributions from iNoutsos et al.l ( 1201 3h and 
iPopov (fe Tu rolla ( 2012) — the emerging magnetic field — is also briefly discussed. 
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1 INTRODUCTION 

Initial parameters of neutron stars (NSs) can be neither observed 
directly, nor calculated from the first principles, yet. Such param- 
eters as initial magnetic field or initial spin period distributions, as 
well as initial velocity distribution, etc. are derived from observa- 
tional data using dif ferent assumptions, somet imes via population 
synthesis mod eling I Popov & Prokhorovll2007l) . 

Recentlv. INoutsos et alj 1 2013 ) (hereafter N13) presented a 
detailed analysis of a sample of radio pulsars (PSRs) with relatively 
well-determined kinematic ages. Assuming standard magneto- 
dipole formula they reconstruct the initial spin period distribu- 
tions for PSRs, and compare it with previous results. In partic- 
ular, thev provide comparison with the distribution proposed by 
IPopov & T urolla (2012) (hereafter PT12) for the population of NSs 
associated with supernova remnants (SNRs) with known ages. 

The sample of pulsars used by N13 has average ages of about 
a few million years. The NSs used by PT12 are much younger The 
two initial spin period distributions do not coincide: the one by N13 
is wider, even PSRs with initial spin periods Po ^ 1 s are present. In 
PT12 the authors provide for many sources just upper limits on Pq 
and they do not provide the initial spin period distribution. So it is 
difficult to compare directly the two distributions. However, as the 
data in PT12 are shown to be compatible with a gaussian distribu- 
tion with (Po) = 0.1 s and crp = 0.1 s, we compare data from N13 
with two gaussians. The Kolmogorov-Smimov test demonstrates 
that the probabilty that the data from N13 are compatible with the 
gaussian with <Po> = 0.1 s and crp = 0.1 is 4.5 lO"'. Slighlty wider 
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gaussians are also in correspondence with the data set in PT12, so 
we also compare the data from N13 with the gaussian (Pg) = 0.2 s 
and <Tp = 0.2 s. The probability that they come from the same dis- 
tribution is 0.0582. We conclude that distributions from N13 and 
PT12 cannot correspond to the same population of sources. 

In this note we discuss how two initial spin period distri- 
butions can be brought in correspondence with each other We 
propose two possibilities: field decay or emerging magnetic field 
buried earlier by the fall-back accretion. 



2 DECAYING MAGNETIC FIELD 

In the paper by N13 initial periods are reconstructed from the 
present ones using the magneto-dipole formula with all parameters 
(magnetic field, angle between spin and magnetic axis) constant in 
time. In PT12 the authors used the same assumptions, but objects 
in their sample are about two orders of magnitude younger For this 
young sample of objects in SNRs, evolution of the field or of the 
angle cannot be very significant, but it can influence the other one 
studied in N13. 

Below we do not speak separately about the angle evolution. 
Instead, one can assume that the "magnetic 'field" in the magneto- 
dipole formula is actually an "effective field", = B sin^ where 
X is the angle between spin and magnetic axes and B is the magnetic 
field at the NS equator 

The process of field decay in pulsars can explain the long ini- 
tial spin periods, inferred under the assumption of a constant mag- 
netic field, as earlier braking was faster than it is now. This would 
widen the reconstructed (under assumption of constant parameters) 
Po distribution. Let us give an example. 

Rapid exponential field decay on the time scale of about few 
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Figure 1. Initial spin period distributions. The liistogramm represents ini- 
tial periods from N13. The solid line is a gaussian with (Pq) = 0.2 s and 
D"p = 0.2 s. The dotted line is a (Pq) = 0. 1 s and crp = 0. 1 s. All gaussians 
are normalized to the same number of objects as in the histogram. 



million years is assumed to be excluded for NSs with standard field 
values (magnetars are not discussed here) basins on the pop ulation 
synthesis of radio pulsars jpaucher-Giguere & Kasp i' 2006). How- 
ever; longer time scales, closer to 10^ years are possible (Pons et aL 
l2009h . In particular, we choose for an illustration a characteris- 
tic time scale for t he exponential field d ecay equal to 5 Myrs, as 
it was proposed by Gonthier et al.l (l2002h (see, however, critics in 
iFaucher-Giguere & Kaspil2006h . More complicated models of field 
decay can also be used. For example, as the NS crust cools down 
a t age s > 1 Myr the Hall drift may become important dPons et al] 



1 20091) : then the rate of decay depends on details of the NS param- 
eters, but for illustration purposes we prefer to use a monotonic 
exponential decay which proceeds with the same rate for all NSs. 

Note, also, that the distribution with (Pg) = 0. 1 s and crp = 
0.1 s, for which N13 make comparison, was given in FT 12 just 
as an illustation, not as the best fit. Variants with larger (Po) 
and (Tp are also possible, as it was noted in that paper. We take 
the following initial distributions. Spin periods have gaussian dis- 
tribution with (Tp = 0.15 and mean value 0.3 s. Initial mag- 
netic field has gaussian distribu tion in log with iTb = 0.55 and 
mean value log £ o/[G]= 12.65 jFaucher-Giguere & Kaspil l2006l : 
IPopov et alj2010l) . 

The magneto-rotational evolution of PSRs in the population 
synthesis model was calculated for exponentially decaying mag- 
netic field: 
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As we want to compare our calculations with results from N13 
then we need to take relatively young (with respect to the whole 
population of normal PSRs) and not very faraway objects. To do 
so we consider only not-too-old pulsars at distances smaller than 
5 kpc from the Sun. The latter condition is related to the potential 
difficulty in estimating the kinematic ages for distant PSRs. In the 
sample from N13 the majority of objects with estimated Pg have 
distances < 5 kpc. In addition, among distant objects dim PSRs 
can be missed, which potentially may result in a bias in the Pq 
distribution. To avoid this, we exclude distant NSs in our model 
limiting the distance to a value close to the maximum in the sample 
from N13. To select young objects we take PSRs with true ages 
below 10' year. On the other hand, since for very young NSs it is 



difficult to estimate kinematic ages we neglect objects with /true < 
10^ yrs (indeed, such NSs are absent in N13 ). 

Let us briefly describe the population synthesis model we use. 
Pulsars are bom with constant rate (the exact value is not impor- 
tant for our study) in four spiral arms which are parametrized by a 
logarithmic spiral (IValleell2005l) . We follow evolution of each pul- 
sar in the model for 5 10^ yrs. The motion of the pulsar is numeri- 
cally integrated in the gravitational pote ntial of the Galaxy. The po- 
tential is chosen in the same form as in lFaucher-Giguere & Kaspil 
( l2006l) . i.e. (f>a(r,z) = (l>dh(r,z) + (pb(r,z) + ^n(r,z)- Here 0dh is the 
potential of the disk and halo, 0b('', z) is the potential of the bulge, 
and (l>„ir, z) is the potentia l of the Galactic nucleus (for details see 
iKuiiken & Gilmorelll989l) . The kick velocity is sampled from the 
double-sided exponential distribution: 
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with (v;) = 180 kms"' l lFaucher-Giguere & Kasp i^'2006^■ and |v/| is 
the absolute value of v/. Every component of the kick velocity is 
randomly generated according to the probability function (Eq.|2ll. 
We neglect the Shklovskii effect as well as changes of the relative 
position of the Sun and spiral arms. 

The dispersion me asure is calculated acco rding to the NE2001 
electron density model dCordes & Lazioll2002l) . The spin period of 
a pulsar at the moment of observation is calculated by integrating 
the standard magneto-dipole formula: 
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Here K = %n^R^I(3Ic^) = lO^^' cm s' g-', is a NS radius, / is 
the moment of inertia, and c is the speed of light. B(t) is calculated 
according to Eq. ([T). 

A pulsar is considered to be detected if its luminosity exceeds 
Sniin for the Parkes multibeam (M anchester et al. 200i)) or Swin- 
burne survey dEdwards et al.l200ll) . and if it is located within a 15° 
(full width) stripe along the Galactic plane. To calculate beaming 
and radi o luminosity we use the same as sumptions as in the best 
model in lFaucher-Giguere & Kaspil ( 12006). A pulsar is assumed to 
be observable till it crosse s the death-line ( Ruderman & Sutherlaridl 
ll975l : lRawlevetal .1119861) : 
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After we obtained a synthetic population of observed radio 
pulsars we reconstruct the initial spin period distribution assuming 
the magneto-dipole formula with constant parameters, i.e. we fol- 
low the procedure used by N13 : 



(5) 



Here P is the present spin period, r is the characteristic (spin-down) 
age, and /true is the true age of a pulsar, as computed from the model, 
in contrast to the approach by N13 where the most probable esti- 
mate based on the kinematic age is used. However, this assumption 
cannot cause problems because typically for PSRs from N13 the 
kinematic age, f^n, should be very close to the true age. 

We use the classical definition of the characteristic age: 
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It is expected that Pq calculated with Eq. (|5) is longer in the case of 
magnetic field decay. Indeed, we can write: 



Pi + 2KB: 
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Figure 2. The histograms show the reconstracted distribution of Pq for the population synthesis model. The population evolved with exponentially decaying 
magnetic field, T,nag = 5 Myrs. Only PSRs with d < 5 kpc (at the moment of observation) are taken into account. The soUd Une in both panels represents the 
true initial distribution of periods used in our calculations. In the left panel we show results where PSRs with characteristic ages > 10^ yrs are neglected. In 
the right panel — objects with true ages > 10^ yrs are not taken into account. In both plots objects with < 10^ yrs are not used. 



where is the initial magnetic field, and f(t) is the magnetic field 
decay function. Then: 



(8) 



where Tcon- = Pl/(2.KBg) is a correction term, and physically this is 
the "initial spin-down age". 

It is known that for normal pulsars (B ~ 10'^ G) with ages 
larger than ~ 10' years Tcon is usually much smaller than r and, 
therefore, changes in the reconstructed Pg in Eq. (|5} are determined 
mainly by the magnetic field decay. In fact, when decay starts to be 
important the second term in Eq. (8) is always larger than the first 
one, and they grow with approximately the same rate. If the field 
evolution, contained in f(t), is the same for all PSRs, then when 
TcoiT starts to be insignificant, the char acteristic ages of PSRs with 
the same ttYuc are approximately equal ( Ilgoshevll2012l) . However, P 
is still determined by So, see Eq. Q. So, the reconstmcted initial 
periods appear to be dependent on Bq. 

Results of population synthesis simulations are shown in Fig. 
[2] As a histogram we present the reconstructed distribution of Pg (to 
be compared with Fig. 17 in N13, see also Fig[TJ. The real initial 
spin period distribution used in the model is shown with a solid 
line. We show also a variant of this plot where the condition /„ue < 
10^ yrs was replaced by t < 10' yrs. The motivation is that in 
N13 all but one PSRs with estimated initial spin periods satisfy this 
condition. As it can be seen, there is virtually no difference between 
two panels. 

Obviously, the reconstmcted distribution looks much different 
from the real underlying initial distribution, but very similar to the 
one derived by N13. This demonstrates that magnetic field decay 
can potentially explain the difference between two distributions for 
realistic parameters. This is the main result of our paper Of course, 
some fine tuning is possible, but we think that with the current low 
statistics this is premature. 



3 DISCUSSION 

3.1 Re-emerging magnetic field 

Another posibility to explain the difference between distributions 
from N13 and PT12 is related to the idea of emerging mag- 
netic fi eld after it has been bu r ied due to strong fall-back ac- 
cretion (lMuslimov&Pagelll995l; IHoI|2011| : IVigano & Ponil2012l : 



iBemal et al.ll2012t 



The initial spin period distribution derived by N13 for val- 
ues below ~ 0.3 - 0.4 s is not much in contradiction with PT12. 
The main problem is due to objects with Pq ~ 1 s, which are just 
a few. One way to explain them is to find a physical reason why 
they are hidden among young sources (analyzed by PT12) and vis- 
ible among the older population (studied in N13). Fall-back which 
buries the magnetic field can do the job if the field emerges faster 
than a few hundred thousand years. 

Among young NSs in supernova remnants there is a group 
of so-called central compact objects (CCOs) for which spin pe- 
riods are not measured as no pulsations (in radio or/and in X- 
rays) are obser ved. They are assum ed to be relatives of so-called 
anti-magnetars dOotthelf eta"i]|2013h which have long Pq: ~ OA - 
0.5 s. Then we can suspect that their initial spin periods are also 
long. One can speculate that they are on average longer than those 
of CCOs with observed pulsations. For example, it is possible 
to propose a hypothesis that for longer initial spin periods fall- 
back is more significant, and so the field is buried deeper, then 
we do not see pulsations, which are due to non-isotropic temper- 
ature distribution produced by the magnetic field. For example, 
iGunevdas & Eksil ( 120131) provide arguments in favour of a corre- 
lation between kick and amount of matter accreted durin g a fall- 
back episode. On the other hand, ISpruit & PhinnevI ( Il998l) discuss 
a correlation between the value of kick velocity and initial spir0. 
Additional kick, in their model, can spin-up the star, so the larger 
the kick, the faster the spin. Then NSs with smaller initial kicks can 



' About fall-back see, for example. fColpi et al .Il996l and references to early 
studies therein. 

^ Note, however, that our analysis using Pq data from N13, PT12, and ve- 
locities from the ATNF catalogue does not support any correlation between 
Pq and velocity. 
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Figure 3. Pq - B plot for sources from N13. Initial spin periods are re- 
constructed from the present day values using the magneto-dipole formula, 
constant field, and kinematic age as the trae age. Aitows point to initial 
parameters of pulsars with reconstracted Po > 0.6 s if the exponential mag- 
netic field decay with Tniag was operating. If for a PSR several arrows are 
plotted, then solid lines correspond to the most probable estimate of kine- 
matic age, from N13, and dashed lines, to the lower and upper limits 
on (kin- If just one solid arrow is plotted, then it coiTesponds to the lower 
limit on because large values do not allow to derive any estimate for 
the initial spin period. 



region in the right panel is undeipopulated as sources are shifted 
toward the bottom-right part of the plot. However, if we remove ob- 
jects with T > 10' yrs (which are absent in the sample from NI3), 
then the plot for the decayed field starts to look relatively similar 
to Fig. [3] If in the original data from N13 the correlation is valid 
for the whole range, then this can be considered as an argument 
in favour of the decaying field model. We additionally illustrate it 
adding to Fig.[3]arrows connecting reconstructed in N13 initial pa- 
rameters with those obtained in the model with decaying magnetic 
field. Clearly, PSRs are moved to much shorter initial periods, and 
so no contradiction with the distribution proposed in PT12 exists. 

In the opposite case, if the correlation in Fig.|3]is just due to six 
objects on the right, then we have a separate population of sources 
with long initial spin periods, and this, in our opinion, can be an ar- 
gument in favour of emerging magnetic field. In PT12 the authors 
did not find any correlation between Pq and B, but in their sample 
of young objects, most probably, there are no (or very few) PSRs 
with re-emerged fields, and magnetic fields could not change signif- 
icantly dueJo_decaY^ote_also, that CCOs can have large toroidal 
fields l lShabaltas &'Laill20I2h . and dipole fields can be at least cor- 
related with toroidal, then objects with long initial spin periods can 
have larger magnetic fields (after they fully emerge), as we see in 
Fig. [3] 

To distinguish between the two possibilities more sources with 
known ages are necessary. 



have longer spin periods and larger amounts of accreted matter. If 
this is the case, then among young NSs (analyzed by PT12 ) sources 
with large initial spins should not be visible, but on the timescale 
10"* - 10^ years their fields can re-emerge (otherwise, they should 
be visible among high-mas s X-ray binaries, but this is not the case, 
see IPopov & Turoiiall2013h . and so they contribute to the sample 
analyzed by N13 . 

3.2 Pq - B correlation 

Here we want to discuss a possible correlation between present 
magnetic field (determined as B ~ VPP) and initial spin period 
reconstructed with Eq. (5). We use the data from Table 1 in N13. 
Results are shown in Fig. |3] It seems, mostly due to a group of 
sources with Po > 0.6 s, that for longer initial spin periods the mag- 
netic field is higher. The plot Po- B presented by PT12 is much dif- 
ferent. So, again we have to explain the differencies. The question 
is: is it a real correlation, or this is an artefact due the assumptions 
made to derive Po, or is it just a fluctuation? 

Statistics is rather poor, and two difl'erent conclusions can be 
made. Either the coiTelation appears only due to the group of six 
PSRs with the longest Pq, or the coixelation is valid for the whole 
range of Po (the two outlying PSRs with the largest fields can be 
due to a fluctuation). Also, of course, it is probable that this feature 
of the distribution is just due to some unknown selecion effects, but 
we do not discuss this possibility further. 

We calculate similar plots with a population synthesis model. 
Results are shown in Fig. |4] In the left panel we show the case 
for the constant magnetic field, and in the right that of a decay- 
ing field with T,„ag = 5 10* yrs. Initial spin periods are recon- 
structed using Eq. ([5}. We select objects with true ages in the in- 
terval 10' < /true < 10' years. Both plots are not similar to Fig.|3] 
no obvious correlation is visible in any panel. The main difference 
between the two distributions is that for decayed fields we have 
objects with long reconstructed Po, as expected, and the top-left 



3.3 Phenomenological magnetic field decay model 

More complicated magnetic field decay models exist, and they also 
can be applied to reproduce the Po - 6 plot with the data from N13. 
As an illustrati on we use a mo del developed by Igoshev, in prepa- 
ration (see also Ilgoshevll20 1 2[) for a different study, but applicable 
also in this case. Parameters and evolutionary laws in the model 
were fitted using the data on r distributions for a large sample of 
known PSRs. 

As there are no pulsars with t > 10' years in the sample from 
N13 it may indicate that (despite the field decay) for ages /n-ue ~ 
10^-10' yrs we have f tme ~ t. The phenomelogical model of decay 
suggested bv llgoshevi j2012l) satisfies this condition. The model is 
parametrised by the following equation: 







a — 




\[ to. 





Parameters of the fine-tuned model have the following values: to = 
10000 yrs, a = 0.034, y = 1.17, C = 0.84. Initial distributions are 
the following: (Pq) = 0.2 s, crp = 0.15 s, and <logSo/[G]) = 12.92, 
cTfl = 0.47. For f(,.„e > 3.5 10' yrs magnetic fields are assumed to be 
constant. 

The distribution of pulsars in Pq - S plane calculated within 
this model is similar to the one from N13, see Fig. |5] Despite 
the fact that this distribution is based on a multi-parametric, phe- 
nomenological model, this exercise illustrates that the field decay 
can be fine-tuned to produce the Pq - B similar to the results by 
N13. 



4 CONCLUSIONS 

In this work we studied if results on the initial spin period distribu- 
tions obtained in PT12 and N13 can be brought in correspondence 
with each other. In both papers the authors reconstructed initial pe- 
riods using independent estimates of NSs ages (SNR age in PT12, 
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Figure A. Pq- B plots for synthetic populations. Initial spin periods are reconstructed according to Eq. (5). Fields represent present day values. Only objects 
with lO' < itfue < lO' yrs ai'e plotted. Left: constant magnetic field. Right: decaying field with T,^ag = 5 Myrs. The solid line in each plot coiTesponds to 
T = lO' yrs, the dashed line in the bottom, to the death line. 
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Figure 5. Pq - B plot for the phenomenological magnetic field evolution 
model. Initial spin periods are reconstructed according to Eq. (5). Present 
day values of the magnetic field are shown. Only objects with 10^ < ?i,uc < 
10^ yrs are plotted. Magnetic field evolution is calculated according to Eq. 
(9). The solid line corresponds to r = lO' yrs. The dashed line in the bottom 
corresponds to the death line. 



and kinematic age in N13) and the standard magneto-dipole for- 
mula with braking index n = 3. One sample (PT12) contains ob- 
jects with ages from a few thousand years up to ~ 100 kyrs. The 
other sample (N13) contains much older NSs with an average age 
of about a few million years. We analyzed if magnetic field evo- 
lution can produce both reconstructed initial period distributions 
starting with the same true distribution. 

We apply the population synthesis technique to scrutinize the 
effects of magnetic field decay on the reconstructed Pq distribu- 
tion. We found that a simple exponential decay with r^ag = 5 Myrs 
can reproduce a narrow distribution for younger objects, and much 
wider one for objects with ages ~ few Myrs, in correspondence 
with PT12 and N13. Potentially, magnetic field evolution can be 
fine-tuned to reproduced additional features of the sample pre- 
sented in N13, but such calculations are beyond the aims of this 
note. 

In addition, we briefly discussed the possibility that re- 
emerging magnetic field can be also used to explain the differences 



between the two initial spin period distributions, as NSs with longer 
initial spin periods can be hidden among the younger population 
studied in PT12. However, we do not model this mechanism. 

We conclude noticing that the differences between initial spin 
period distributions obtained by PT12 and N13 are due to evolution 
of the magnetic field on the time scale smaller than a few million 
years. To distinguish between different models of field evolution it 
is necessary to increase statistics of NSs with known ages. 
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